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Abstract 

We discuss the solitary wave solutions of a particular two-component scalar field model in two- 
dimensional Minkowski space. These solitary waves involve one, two or four lumps of energy. The 
adiabatic motion of these composite non-linear non-dispersive waves points to variations in shape. 

1 Introduction 

A solitary wave travels "without changing its shape, size, or, speed' , pQ. In this paper we shall deal with 
solitary waves in relativistic scalar field theory in a two-dimensional space-time. Besides the prototypes, 
the kink of the non-linear Klein-Gordon equation and the sine-Gordon soliton found in models where 
the scalar field has only one component, solitary waves have also been discovered in systems with two- 
component scalar fields, j2]-jSl- 

In a deformed 0(2) Linear Sigma model, the celebrated Montonen-Sarker-Trullinger-Bishop (MSTB) 
model, see [I], [S], there exist solitary waves with two non-null field components. The model was first 
proposed by Montonen 0], and Sarker, Trullinger and Bishop (Sj in the search for non-top ological solitons. 
Nevertheless, in both papers the existence of two kinds of solitary waves, respectively with one and two 
non-null components, was noticed. Slightly later, Rajaraman and Weinberg [2] proposed the trial-orbit 
method to study the two-dimensional mechanical problem equivalent to the search for static solutions. 
By these means, they found that the TK1 -one-component topological- kink is given by a straight line 
trajectory in field space, whereas the TK2 -two-component topological- kinks come from semi-elliptic 
trajectories; they also enlarged the list by discovering NTK2 -two-component non-top ological- kinks 
bound to elliptic orbits. Through numerical analysis, [H], [I], Subbaswamy and Trullinger, looking for 
more exotic orbits, observed the existence of a one-parametric family of NTK2, including that previously 
discovered in [2]. Moreover, they discovered an unexpected fact; the "kink energy sum rule": NTK2 
energy is the same as the addition of TK1 and TK2 energies. Magyari and Thomas, [5J, realized that 
the reason for the sum rule lies in the integrability of the equivalent mechanical problem. Ito went 
further to show that the mechanical problem giving the solitary waves as separatrix trajectories in the 
MSTB model is not only integrable but also Hamilton- Jacobi separable, and, moreover, distinguished 
between stable and unstable solitary waves, see by applying the Morse index theorem. The full 
Morse Theory for the MSTB model was developed by one of us in 

Only TK2 kinks are stable and genuine solitary waves in the sense that they travel without distortion 
in shape, size, and, velocity. In a model proposed by Bazeia-Nascimento-Ribeiro and Toledo -henceforth 
the BNRT model, |12j . \1'A\ - things are different. In this case, both the TK1 and the TK2 kinks -discovered 
in the papers quoted above- are stable and degenerate in energy. In fact, they are distinguished members 
of a one-parametric family of kinks found in |14j . all of them degenerate in energy, but composite in a 
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certain sense. Apart from the center of the kink, there is a second integration constant; for large values 
of the second integration constant b solitary waves seem to be composed of two lumps of energy, whereas 
for b small they only show one lump. In References JT5] and Pj5] we first unveiled the whole manifold 
of solitary waves of this model and then showed that in this system solitary waves can travel without 
dispersion, although changing the relative position of the two lumps. 

The distribution of the lumps is, however, completely symmetric with respect to the center of mass. 
In this paper, we shall study a similar system to the MSTB and BNRT models in this framework. The 
model can be understood as the dimensional reduction of the planar Chern-Simons-Higgs model, |19j . to 
a line and zero gauge field. The whole manifold of composite solitary waves was described in Reference 
|17j . We shall again derive the solitary wave solutions of the model called A in the Reference cited, now 
using the Bogomol'ny method, see JH]- This strategy is better suited for describing the wave properties 
than the solution of the Hamilton-Jacobi equation. Later, Manton's adiabatic principle, designed to 
elucidate the slow-motion dynamics of topological defects, will be applied to our solitary waves. 

The organization of this paper is as follows: In Section §.2 we introduce the model and describe the 
spontaneous symmetry-breaking scenario. We also identify two systems of first-order differential equations 
that are satisfied by solitary waves using the Bogomol'nyi procedure. In Section §.3, the solitary wave 
solutions of the first-order equations will be obtained. We shall show that there exist two basic kink- 
shaped solitary waves and that the rest of the kinks are configurations of two or four basic kinks. The 
problem of the stability of these waves will be also addressed in this Section; the double solitary waves 
are stable whereas quadruple kinks are unstable. In Section §.4, we apply Manton's principle to study the 
adiabatic dynamics of solitary waves as geodesic motion in their moduli space. There is no problem with 
double solitary waves but we also analyze the motion of quadruple lumps by pushing them in the direction 
of instability in configuration space. Finally, we explore how the system react to a small perturbation of 
the potential term, preserving the basic kinks; the induced forces give rise to a remarkable bound state 
of kinks, similar to the breather mode of sine-Gordon theory. 



2 The model: quintic non-linear Klein-Gordon equation 



We shall focus on a (l+l)-dimensional relativistic scalar field theory model, whose dynamics is governed 
by the action: 
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and xiUfi) = XiiU/i) + ^X2{y^) '■ R — * C is a complex scalar field. We choose goo 
9i2 = 921 = as the metric tensor components in Minkowski space R > , whereas 
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sets the non-linear interactions in A, m and (5 are coupling constants of inverse length dimension. 



Introducing non-dimensional variables x 

■2 



f4>, = and p- = a 2 , the action functional reads: 
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— V((f), <jf) -plotted in Figure la- is a semi-definite negative polynomial expression of degree six which 
depends on the unique classically relevant non-dimensional coupling constant a. 
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Although the potential that we have chosen seems to be very awkward it is physically interesting 
for two reasons. First, V(<pi,<p2) is a deformation of Vcsh(</>) = 2~ (p* $ (4>* 4> ~~ I) 2 that appears in the 
celebrated Chern-Simons-Higgs model, see ^H] and References quoted therein. The CSH model is a (2+1)- 
dimensional gauge field theory for a complex scalar (Higgs) field where the kinetic term for the gauge 
field in the Lagrangian is the first Chern-Simons secondary class, whereas the Higgs self-interaction is 
ruled by Vcsh- A phenomenological theory for the fractionary quantum Hall effect has been established 
on the basis of a non-relativistic version of the CSH model in [20]. Dimensional reduction to (1+1)- 
dimensions at zero gauge field leads to a class of systems including our model because there are no 
Goldstone bosons on the line [21]; the infrared asymptotic behavior of the quantum theory would require 
the modification of Vcsh in such a way that the zeroes become a discrete set and massless particles 
are forbidden. The deformation of our choice comply with this requirement and, moreover, ensures 
integrability of the dynamical system to be solved in the search for kinks. Second, we shall show in 
Section §3 that the Hamilton characteristic function, the "superpotential" W 1 ^) (|16[). is precisely the 
potential of the MSTB model. This fact links both systems in a hierarchal way. 

The second term in (J2J) breaks the U(l) symmetry of this model. The G = Z2 x Z2 group, generated 
by the reflections (pi — ► —(pi and <p2 — ► —<p2 in the internal space C, leaves, however, our system invariant. 




Figure 1: a) Potential term —V((p,(p*). b) Set of zeroes M of — V((p, < 
The field equations form the following system of coupled second-order PDE: 



d 2 cpi d 2 



dt 2 dx 2 
d 2 6 2 d 2 6 2 



i[l + 3(Pf + 34 + (P 2 (Q(P 2 - 4) + 2(a 2 - 2)<t>l] (3) 
2 [3(P\ + 2 ( p 2 {3 ( p 2 + a 2 -2) + ((P 2 + a 2 -\){3(P 2 + a 2 -I)] . (4) 



dt 2 dx 2 

The PDE system ©, Q is akin to the non-linear Klein-Gordon equation with quintic -besides cubic- 
non-linear terms. Initial, (p(to,x), j^(to,x), and/or boundary conditions, 0(t,±xo), ^(t, ±xq), must be 
chosen in order to search for physically relevant solutions. We focus our attention on those solutions that 
can be interpreted as lumps or extended particles: that is, solitary waves or kinks. Recall the definition 
of a solitary wave, see e.g. 0: "A solitary wave is a localized non-singular solution of any non-linear 
field equation whose energy density, as well as being localized, has space-time dependence of the form: 
e(t, x) = e(x — vi), where v is some velocity vector". 



2.1 Structure of the configuration space 

The set of zeroes of V((pi, (P2) -maxima of —V((pi, (P2)- is: 

M = {4> A± = ±1; 4> B± = ±ia; 4>° = 0} , a = \/\-a 2 

see Figure lb. These constant configurations are the static homogeneous solutions of the PDE system 
©, Q because they are critical points of V(<pi, (p2)'- ■§¥- = = J^-. Note that the reflection (pi — > —<f>\ 
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sends <p A+ to (p A and vice- versa; they belong to the same G-equivalence class: <p A = {<p A+ , (p A }. There 
are another two G-orbits between the constant solutions: <ft B = {(p B+ , (p B } and <ft°. The moduli space of 
homogeneous solutions- i.e., the set of zeroes of V modulo the symmetry group- consists of three points in 
C: Ai = =|=r = {<p A , (p B , 4>°}- The degeneracy of the homogeneous solutions -all of them have zero energy- 
causes spontaneous breaking of the discrete symmetry. The constant solution cp A breaks the symmetry 
under G = 7Li x 7L<i transformations of the action functional Q to the little group Hi = {e} x Z2 
generated by <p2 —<p2- Simili modo, the remaining symmetry group at <p B is the H2 = Z2 x {e} little 
group generated by (pi — > —(pi, whereas the point (p preserves the full symmetry G. 

The configuration space of the system C = {(p(t,x) G Maps(M, C) / £[(p] < 00} is the set of maps 
from 1R to C for fixed time t such that the energy functional 1 
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is finite. Thus, every configuration in C must comply with the asymptotic conditions: 

lim (f>(t,x)eM lim M^ = . (6) 

z^ioo x— >±oo OX 

These asymptotic conditions play an important role from a topological point of view. The configuration 
space is the union of 25 topologically disconnected sectors: 

C= C AA U u C°° u C AA u Cg 
u Cil u Cli u C AB u Ci A 

u Ci° u Cg A u C BO u C^ B (7) 
The relation of this notation to the asymptotic conditions © is self-explanatory. Three examples: 

1. Sector C°° . Boundary conditions: 

lim (p(t, x) = 4>° , lim (p(t, x) = 0° 

X— ► — 00 x— >oo 

2. Sector C AB . Boundary conditions: 

lim (pit, x) = 4> A± , lim (p(t, x) = 4> B± 

x— >— 00 x— >oo 

3. Sector C BA . Boundary conditions: 

lim (p(t, x) = 4> B± , lim (p(t, x) = (p A ^ 

X — > — 00 X — *oo 

Because temporal evolution is continuous (a homotopy transformation), the asymptotic conditions © 
do not change with t and the 25 sectors ((7J) are completely disconnected. Physically this means that a 
solution in a sector cannot decay into solutions belonging to other different sectors; it would cost infinite 
energy. 

The five homogeneous solutions of (EI)-® that are zeroes of V belong to the five sectors C±±, C BB 
and C°° respectively. Small fluctuations around one of them, t()(t,x) = cp + 5ip(t,x), (p G A4, solve the 
linear PDE system: 

2 - - d 2 V 

(VS ab + M 2 a S)) SMt, x) = , M 2 ab (4>) = q^JuS^ (8) 
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1 Strictly speaking the energy is A 2^/2 ^1 as gi ven m formula ©, is a non-dimensional quantity. 
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The solution of (JHJ) via the separation of variables -8tp a (t, x) = exp[iu;t]f^(x)- leads to the spectral problem 
for the second order fluctuation -or Hessian- operator: 

i + mis) \ ( m*) \ =t *( fi( x ) \ 

mis) -£+m$s) ) v m*) ) v m*) ) 

Ti.{4>) is a diagonal matrix differential operator for the three points of M with a positive definite spectrum; 
every constant solution belonging to M. is stable. Thus, there are three types of dispersive wave-packet 
solutions, living respectively in one of these sectors. The building blocks are plane waves with their 
dispersion laws respectively determined by: 

M W ) = ( I ° ) M W = ( ; 4 4 » 4 ) M HiP) = ( J » ) . 

The two branches for each type are: 



<j) A 




0° 


uj 2 (k) =k 2 + 4 
u \ q ) = q 2 + a* 


o; 2 (A;) = P + a 4 
w 2( g ) =g 2 +4a 4 


Cj 2 (k) = k 2 + 1 
o; 2 (g) =g 2 +a 4 




2.2 Solitary waves from integrable dynamical systems 



Any other static solution of ©-(JU is a solitary wave that lives in one of the remaining twenty topological 
sectors of the configuration space. The field profiles can be either kink- or bell-shaped, and Lorentz 
invariance allows the use of the center of mass system. The energy density and field equations read: 



d 2 



dV 



dx 2 



1,2 



(9) 



The search for solitary waves is tantamount to the search for finite energy static solutions cf>(x) of the 
PDE system ©-©, which, in turn, reduces to the ODE system ([9*]). A Lorentz transformation sends the 

static solution (j){x) to <f>(t,x) = (j) ( J 1^2 ) • 

The ODE system (jHJ) is nothing more than the equations of motion for a two-dimensional mechanical 
system: understand (^>i,<fo) as the "particle" coordinates; x as the "particle" time, and U = —V as the 
"particle" potential energy. We shall use a mixture of Hamilton-Jacobi and Bogomol'nyi procedures to 
solve this system. 



dW 



dW 



(10) 



is the "time" -independent Hamilton-Jacobi equation for the mechanical system and zero "particle" energy. 
Using the Hamilton characteristic function W, the solution of (|10[). we write the field theory potential 
energy for static configurations a la Bogomol'nyi: 



£[<j>i,<h 



dx 
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dx 
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dx 
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(12) 
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is a topological charge, depending only on the sector of the configuration space if W is well-behaved 
enough to apply Stokes theorem: 



T = W(fa(oo), Moo)) - w(M-oo), M-oo)) 

In this case, an absolute minimum of the field theory energy fil[). in the corresponding topological sector, 
is reached by the solutions of the following system of first-order equations: 

d^i = dw_ dfo = dw 

dx d(pi dx d(p2 

(in the squared terms are always positive and the other is a topological constant). It is easy to check 
that solutions of the first-order equations (fT3)) also solve the second-order (Euler-Lagrange) equations © 
and consequently © and (@J). The energy of these solutions depends only on the topological sector of 
the solitary wave, and it saturates the topological bound £[</>k] = The quotient of the first by the 
second first-order equations gives the "orbits": 



dW , , dW 



(14) 



( dW (dW ~ )^ 

<,-QT^{fa,fa(fa))> dfa = dx , | — (fa(fa), fa) J d(p 2 = dx , (15) 

where fa(fa) (respectively fa(fa)) is a solution of (|14|). provide the "time" -schedules of the "particle". 
In fact, this system of three equations would be also obtained in the Hamilton-Jacobi framework setting 
all the separation constants equal to zero. 



3 The variety of solitary waves 

Using elliptic coordinates, 

u{x) = ~ (j \fa{x) + a) 2 + (f)l{x) + ^/(0!(x)-a) 2 + 02(x)^ E (a, +00) 

v(x) = l - U \m%) + °Y + 4>l{x) ~ V '(M*) ~ v? + G (-^ a) 

one can show that the equation ()10j) is separable for our choice of V(4>\,(p2)- Thus, two independent 
solutions for W as a function of u, v are found. Back in Cartesian coordinates, these are: 

W I (<t>) = ^[(4> 2 1 + 4-l) 2 + 2a 2 4] , (16) 

WH W = \ \/(fa + + & yj{fa ~ v) 2 + ~k (<f>l + 4 + ° 2 -2) ■ (17) 

W 1 is a polynomial function of fa, fa without singularities and everything proceeds smoothly. The 
partial derivatives of W n , however, are undefined at the points (cj)\ = ±cr, 4>2 = 0); it is uncertain how to 
interpret this second solution of the Hamilton-Jacobi equation. Note also that if a function W(fa verifies 
formula 1)10(1 then —W(4>) also satisfies this condition. Replacing W by —W merely flips the sign of the 
topological charge. 
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The system of first-order differential equations determined by the first solution W 1 is: 



dx 



dx 



(18) 



whereas the flow of the gradient of the second solution W n is ruled by the first-order ODE system: 
#x _ <M</>f + 4 + a 2 + $(2$ - 1 - a 2 ) + <j>l(a 2 - 1)] 



dx 



dx 



M0i + 4 + <fi(Wl - 1) + ^V 2 - 1) + <t>li^ 2 



i)] 



(19) 



v /( 01 _ fT )2 + ^ v / (01+fj) 2 + ^ 



3.1 Solitary waves: basic types 

In order to solve the ODE systems ([18ft and IJ19JI we start by looking for restrictions on <f> that make both 
systems coincide. The idea is akin to the Rajaraman trial-orbit method, see [3|: guess a trajectory that 
joins two points in M. and minimizes E. It happens that there are two types of curves in C for which 
lfTH|) and (|TU)l become identical; they provide the basic solitary waves of the model. 

1. K.f°: We shall now try the condition fy\ = in systems Q18J) and/or (|19|) . We find kink solutions 
joining the points O and i? in the moduli space: 



4>(x;a) 



(7 



(-l) a i-^\/l +(-l)/ 3 tanh[a 2 (x 



a, f3 = 0,1 



(20) 



Because the orbits stay on the straight line v = in the elliptic plane joining B with O, we shall 
refer to these four solutions 2 -here we consider the parameter a to be fixed- as Kf°. A solitary 
wave connects the points O and B + if a = /3 = and the points O and S_ if a = l,/3 = 
-living respectively in C± B -. For the cases a = 0, /3 = 1 and a = /? = 1, we find respectively the 
antikinks of the above mentioned solutions. One of these solutions has been depicted in Figure 2, 
in addition to its density energy and its orbit in the internal complex plane. In Figure 2c, we see 
why these solitary waves are interpreted as lumps of energy or extended particles. The energy of 
these solutions is: 
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Figure 2: Kf° Kink: a) Orbit, b) Kink profiles and c) Energy density. 



2 The subscript refers to the number of lumps, which is equal to the number of pieces of straight lines that form the kink 
orbit in elliptic coordinates, see |17) . 



2. K AB : On the ellipse 



1 -u = 1 in the elliptic plane - (see Figure 3a), equations (|T%|) and 



also coincide. The following eight solitary wave solutions are found: 



<f>(x;a) 



l) c 



V2 



l + tanh[(-l)/ 3 a 2 x] + i 



V2 



1 -tanhf(-l)/ 3 cj 2 x] 



(21) 



where a, (3, 7 = 0, 1 , x = x — a and a is fixed. All of them join the points (f) A and 4> B in Ai and 
thus belong to the C AB sectors. We shall denote these solitary waves depicted in Figure 3b as K AB 
because their orbits lie on the straight line u = 1 . If = 7 = 0, solution (J21[) goes from B + to 
A + when (3 = (see Figure 3b) and from A + to B + when (3=1. The change in the value of a is 
tantamount to a reflection <$>\ — * —<j)\. Thus, for a = 1, 7 = (3 = 0, (|21|) is a solution running from 
B + to i4_. Analogously, swapping the value of 7 in (|21|) is equivalent to a reflection of the <p2 axis. 
For instance, (|21j) . with a = (3 = and 7 = 1, is the solution that starts from !? + and ends in A-. 
Changing 7 kink and antikink are exchanged. 
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Figure 3: K^ B Kink: a) Orbit, b) Kink profiles and c) Energy density. 

The energy density of these solitary waves is localized (see Figure 3c) in a single lump (thus, K^ 5 ) 
Their energy is: 



S[K AB ) 



\W 
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3.2 Double solitary waves 

There exist many more solitary waves in this model than those described in the previous sub-Section. 

For example, we can try the straight line on the 4>2 = axis joining the points A and O of the moduli 
space Ai, see Figure lb. Plugging this condition into both ODE systems, (fT%|) and ((T§|) . we find the 
common solution: 

4>{x;a) = {-l) a ^=Jl + (-l)^tanh(x - a) a, (3 = 0,1 . (22) 
a/2 v 

In fact, setting the solitary wave center at a, we have found four solutions: (1) a = (3 = 0, the trajectory 
runs from the point O to A + as x goes from — oo to oo; it belongs to C+ A . (2) a = 1, (3 = 0, the trajectory 
starts from the point O and arrives at A^, living in C® A . (3) a = 0,(3 = 1 provides the solution in C AO . 
(4) a = 1,(3 = 1 leads to C A °. The energy of these solutions is: 

£[K£ ] = \T\ = \W I [t A ]-W I tt ]\ = ± . 

From the profile (1221) of the real component of the field we see that these solitary waves are kink-shaped, 
interpolating between the homogeneous solutions (j) A and (fP . In Figures 4(center) and 5(center) we 
respectively depict this solution and its energy density. 

In the usual convention for kinks these solutions are of the type TKl AO , only one-Cartesian component 
being different from zero. However, we denote these kinks as K AO because they consist of two lumps, or, 
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equivalently, their associated orbits in the elliptic plane are formed by two straight lines: (a) from A to 
F along v = ±<r; (b) from F to O along u = a, where F stands for the common foci of the ellipses and 
hyperbolas defining the orthogonal system of coordinates, see |17j . 

There are many more solutions in the same topological sectors as (|22[), In |17j it is shown that 
system (|18[) also separates into two independent first-order ODE when written in elliptic coordinates. 
The general solution can be found and translated back to Cartesian coordinates. The solutions "live" on 
the orbits: 

a^l = b\amf[a 2 (l-4>l)-4>lY , (23) 

which are solutions of equation (fl4*|) for W = W 1 , see Figure 6(left). The explicit dependence on x is 
given by solving (fT51) : 



(x;a,b) 



+ i 



ba 



+ a 2 e 2(x-a) + h 2 e 2^(x-a) ^ + a 2 e 2^{x-a) + & 2 e -2^{x-a) 



(24) 



where a, b £ R are two real integration constants. Constant a sets the center of the solitary wave. The 
meaning of b is two- fold. On one hand, b parametrizes the family of orbits (j23j) linking <f> A and <f>° , see 
Figure 6(left). On the other hand, b enters the solitary wave profile ()24j) . shown in Figure 4 for a = 
and several values of b. For high values of |6|, the real component of the field has the shape of a kink 
but the imaginary component is bell-shaped. When \b\ diminishes, the imaginary component diminishes 
until it reaches zero for 6 = 0. A parallel plot of the energy density for the same values of b shows two 




L 




Figure 4: Members of the family Kg (b): from left to right b = —45, —2, 0, 2, 45. 
energy lumps for high values of |6| and only one lump for |6| close to zero, see Figure 5. We denote these 



Figure 5: Energy density of the previous kinks: from left to right b = —45, —2, 0, 2, 45. 

solitary waves as K^°(b) and realize that the K^°(b) kinks seem to be made from the two basic kinks 
K.f° and K^ B whereas K^°=K^°(0). Both facts are confirmed by the following tests: (1) the \b\ = oo 

limit in (|23|) is only compatible with (f>\ + % = 1 and (pi = whereas the 6 = limit provides us with 
the condition = 0. (2) There is a new kink energy sum rule for all 6: 



AO, 



\W 



W\4P)\ = £[K^ U ] + £\Ki Aa ] = - A 



BOi 



ABi 



1 



(25) 



As a curiosity, for the values of the coupling constant o~ 2 = \ and 6-parameter 6 = ±1 the K^°(±l) kinks 

= (i) 2 .The 



live on semi-circles of radius \ with their centers on the points (±^,0), i.e., (i^i ± ^) 2 + - 
profile is: 



ix;a) 



±i ( 1 +tanh^i) ±2isech , 



(x— a) 
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Figure 6: (left) Orbits for Kf° with 6 = 2 (dotted line), b — 1 (dash-dotted line), b = 0.5 (broken line) and b = 
(solid line), (right) Kj 35 kinks with b — 7 (dotted line), 6 = 3 (broken line) and 6 = (solid line). 



3.3 Quadruple solitary waves 

In Reference |17| we solved again profiting from the fact that this coupled ODE system becomes two 
ordinary uncoupled differential equations in elliptic variables. Thus, it is possible to obtain the equations 
of the orbits and time-schedules explicitly. Translating the results back to Cartesian coordinates, we 
found the complicated expressions: 



~-b 2 [a 2 -C$i$ 2 ) + D(fa,fa 
[a 2 -C{$ u fa) + D(fa,fa 
b 2 [C$i, -a 2 



(C$i$ 2 ) - D$ lt 4> 2 )f - C$i, h) + D{fa,h) 
(C$i, h) + D{4> x ,h)f - C$i, h) - Dtfufa) 
(C$i, h) + D(h,h)f - C$i, fa) - £>$i$ 2 ) 
(C$i$ 2 )- J D$i,^ 2 )f 2 -C$ 1 ,^ 2 ) + £>$i,0 2 ) 



(26) 
,(27) 



where 



C$ 1 $ 2 ) = i$? + ^| + a 2 ) 



1 



$2 + ^ )2 + 2a2( ^2_^2 )2 + (T 4 



as solutions of equation (|TH) for = see Figure 6 (right). 
Also, we showed that 

(-l) Q Cl(l + e 2^ 2 ((^)+^ 2 )) 



4>i(x;a,b) 



4> 2 {x;a,b) 



y/ a 2 ZjZ ^2 e 2(a;-a) + e 2a 2 ((s-a)+bff 2 ) y 7 ! ZjZ -2g2CT 2 ((x-a)+fe CT 2 ) + & 2 e ~2a 2 ((x-a)-ba 2 ) 



aa 2 e ba (e 2( - x - a 1 - 1) 



yV 2 _|_ ^2 e 2(a;-a) _|_ e 2a 2 {(x-a)+ba 2 ) e 2a 2 ({x-a)+ba 2 ) _|_ a 2 e 2((x-a)+ba 2 ) _|_ 



2 26cr 2 



(28) 



solve the second-order equations of motion. 

Besides the continuous parameters a, b £ M, the family of solitary waves (|28|) depends on the discrete 
value of a = 0,1. All the solitary waves in the family ()28|) belong to the sector CE+', their associated 
orbits go from (j) B ~ to B+ when x — a goes from —cxd to +oo. When x = a the foci (j> F± = ±<r of the 
ellipse are respectively reached for a = and a = 1 3 . The foci are the zeroes of -D$i,^ 2 ); moreover, 
the half-orbit equations ((2*6^1 and (|2Tj) are satisfied at these points by all values of b, = b 2 ■ 0. Therefore, 
an infinite number of these trajectories arrives in and leaves from the foci of the ellipse at the "instant" 
x = a, see Figure 6 (right) and note that 

#2 _ + $? + 4>l- 1)$? + 4>l) + ^ 2 (20l - 1)) 



h Mv 2 (-$ + $ + + (0? + $ - 1)$? + 0!)) 



^Solitary waves of the family 1281 with a = always remain in the half-plane cf>x > 0. If q = 1, the solitary waves live in 
the half-plane 4>\ < 0. 
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is indeterminate at <b F± = ±cr, although taking the limit we obtain ~~ = ± 2eT ° at the foci. 

Related to this fact, we stress the following subtle point: the solitary wave family ()28|) does not 
strictly solve the first-order ODE system ()19() because the associated orbits satisfy (|26|) if x £ (—00, a) 
and (|27f) if x £ (a, +00). Taking the 6 2 — ► 00 limit in ()26J) . we find the four K^ fi kinks and two times the 
AF pieces of the orbits; all the orbits in (|26|) have cuspidal points at the foci. In the b 2 — > 00 limit 
of (J2ZJ), the two Kf° kinks are obtained twice, together with the FO pieces of the orbits; curves in 
this family also have cuspidal points at the foci. Smooth trajectories require the glue of half the orbits 
in (|26|) to their complementary halves in (|27|) . This means that these solitary waves are solutions of ()19() 
for —W 11 if x £ (—00, o) and properly for W 11 if x € (a, +00). 

Nevertheless, they are bona fide solitary wave solutions of the second-order field equations although 
the Stokes theorem apply only piece- wise to the integration of dW 11 along these orbits. The energy of 
solitary waves of this type is not a topological quantity: 

£{Kf B {b)) = \W I \^ B -)-W I \<p F± )\ + \W I \<t) F± )-W I \^ B +)\ = ^ . (29) 

Because of the energy sum rule, 

S{Kf B {b)) = £{K BO ) + £{K$ A ) + £(Kf B ) = 2£{K^ B {b)) + 2£{K BO ) , (30) 

we call these solitary waves, composed of four basic kinks, K BB . Replacing x by — x in (J28|). similar 
solitary waves are found -K BB antikinks- living in C BB . 




Figure 7: Members of the family K B : from left to right b = —20, —7, 0, 7, 20. 



e U> °' L ' £U) "^J e(a1 ' J six) " eu) 




Figure 8: Energy density of the above kinks: from left to right b = —20, —7, 0, 7, 20. 

In Figures 7 and 8 we respectively plot the kink profiles (|28ft and energy density for the b = —20, 
—7, 0, 7, 20 members of this solitary wave family. The main novelty is that the real component has the 
shape of a double kink, a fact clearly perceived from the b = ±20 plots; the imaginary component is 
bell-shaped (slightly distorted). The energy density graphics show these solitary waves as configurations 
of three lumps. We easily identify two of these lumps: the lower lump corresponds to the basic K AB kink 
whereas the next lump in height is indeed one of the basic K BO kink, compare Figures 3 and 8. The 
third lump, always in between of the other two, is precisely the K^°(0)= K^ kink. Taking into account 
that the latter are composed of the other two basic kinks, each member of this family is a configuration 
of four basic lumps, two K AB and two K BO kinks, with the peculiarity that one K AB and one K BO are 
always superposed for any value of the parameter b. In particular, the four lumps coincide at the same 
point for the kink configuration characterized by the value 6 = 0. This picture completely explains the 
intriguing kink energy sum rule 1)30(1 . 
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For the a 2 = | and b = ±1 values, the solutions (|28|) live on semi-circular orbits of radius centered 
at the origin of the internal plane C. The analytical expressions are very simple: 



V2(l + 2e( x ~ a )) 



2y/2e^ +i[l-2e {x - a) 



3.4 Stability and gradient flow lines 

The second-order fluctuation -or Hessian- operator 



H[4>(x; a, b)] 



■£+M 2 S(x;a,b)} Ml3(x-aM 
Ml [4>{x; a, b)} + M| 2 $(x; a, b)} 



valued on a solitary wave solution <fi(x; a, b) is a matrix differential operator of Schrddinger type. Negative 
eigenvalues appear in the spectrum of H when a solitary wave is unstable against small fluctuations. 

In general the spectral problem of 7i is non-manageable. There is always, however, important infor- 
mation available: to each one-parametric family of kink solutions 4>(x, c) is attached the eigenfunction 
^ that belongs to the kernel of Ti, i.e., ^ is a zero mode. This is easily checked by taking the partial 
derivative of the static field equations Q with respect to c. 

Stability of Kf° kinks: The Hessian operator valued on the solution (|20l) is: 



BO\ 



Hn ^^/-^-^f^ + M^tanhx + ^tanh 2 , \ 

U ^ W22 / V -in+v* [-f + (-l) /3 |tanhx + ^ tanh 2 x] / 

(31) 

In this case, the Hessian reduces to two ordinary Schrddinger operators of Posch- Teller type. Tin 
rules the orthogonal fluctuations to the solution in the internal plane, whereas H.22 takes into account 
the tangent fluctuations to the K.f° kink. There are no discrete eigenvalues in the spectrum of lri\\. The 
continuous spectrum starts at the threshold cr 4 ; ui 2 (q) = q 2 + cr 4 is non-degenerate in the [<7 4 , 1] range 
but doubly degenerate for eigenvalues higher than 1. H.22, however, presents a discrete eigenstate with 
eigenvalue ui 2 = (zero mode) and eigenfunction ^(x, a). The continuous spectrum uj 2 = q 2 + cr 4 is 
non-degenerate for uj 2 G [a 4 , 4a 4 ] and doubly degenerate for uj 2 > 4<r 4 . The spectrum of 7i{Kf°) is semi- 
definite positive and the Kf° kinks are stable. The zero mode identified is associated with translational 

d£. 



invariance; therefore, the perturbation ip^—Q = j$-(x,a) on the kink causes a infinitesimal translation of 
the lump. 

Stability of the K^ifi) kinks: The Hessian operator valued on the K^°(0) kink is also diagonal: 
*>ffM I H ^ Wl2 A (-£-! + (-l) /3 |tanhx+ f tanh 2 x 



^21 n 2 2 J \ o _j^_ a 2 & 2_ 1 +( _ 1) /3( a 2_ i)tanhx + |tanh 2 x / 

(32) 

Tin, in this case ruling the tangent fluctuations on the solution, presents a zero mode uj 2 = 0, with 
eigenfunction 

^SLo = |^ (^; a, 0) oc - 1 , (33) 

oa J2cosh 3 (x-a)e( x - a ) 

associated with translational invariance. The continuous spectrum shows non-degenerate states in the 
range J 1 G [1,4] and doubly degenerate states for uj 2 > 4. H.22, determining the behaviour of the 
orthogonal fluctuations to the solution, admits another zero mode. The eigenfunction 

4 2 ) = ^r(x;o,0) occosh-h(x-a)e^- a2 >- a) (34) 
ob 
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belongs to the kernel of H22 and exists because the K^°(0) kink is a member of the one-parametric 
family K AO (b). Applying the perturbation (jUJ) to the solution K AO (0) of we obtain K^°(<56). The 
continuous spectrum is non-degenerate in the range 10 2 G [ct 4 ,<t 4 ] and doubly degenerate for > a 4 . 
Therefore, because of the positive semi-definiteness of the spectrum of Ti we conclude that the K AO kinks 
are stable against small perturbations. 

Stability of the K AO {b) kinks : The Hessian for the other members of the K AO (b) family with b 7^ 
is not diagonal. Only the kernel of Ti is known, spanned by the translational zero mode ^(x;a, b) and 

the distortional zero mode ^(x;a,b), orthogonal to each other. Both eigenfunctions obey a neutral 
equilibrium around a given kink configuration. 

It is possible, however, to conclude that the remaining eigenvalues are positive by indirect arguments. 
The solutions to the first-order system are the flow lines induced by the gradient of W 1 . The Hessian 
matrix 

ewa) $£-(</>) 

valued at 4> A , <j) B and <f>° is respectively: 

hess'(^) = (j ° 2 ) , hess'(^) = ( ~f 2 ° 2 ) , hess^°)=(- 1 _°_ 2 

Thus, <f> A is a minimum, cj) B a saddle point and (jP a maximum of W 1 , see Figure 9a. The flow of 
gr&dW 1 runs from (f> A to <jP along the non-intersecting orbits (|23|) . There are no focal/conjugate points 
and Morse Theory ensures that all these kinks are stable, see ^H] and 



d 2 W r rZ\ d 2 W' , 




-0.21 
-1 



Figure 9: a) W z (<h.,<h), * = 75- b) W n (<f>i,<h), <? = ^- 

Stability of K AB kinks: The K AB orbits flow from a maximum (4> A ) to a saddle point (<fi B ) of W 1 . 
They appear at the b 2 — > 00 limit of the K^ family, no dangerous points are crossed, and the Morse 
index theorem tells us that the K AB kinks are stable. 

Instability of the K BB (b ) kinks: Ti valued at the K.f B kinks is again a non-diagonal matrix Schrodinger 
operator; there is analytical information available only for the kernel of Ti. The distortional zero mode 
^(x; a, b), a Jacobi field in the language of Variational Calculus, is zero at the foci of the ellipse for each 

member of the family (|28j): ^(x = a; a, b) = 0. This means that there exists at least a negative eigenvalue 
in the spectrum of Ti because the ground state of (sufficiently well-behaved) Schrodinger operators has 
no nodes. Therefore, the K BB (b) kinks are unstable. A question arises: what is the fate of the K BB 
solutions when a perturbation associated with the negative eigenfunction is exerted? We shall address 
this topic later. 
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A better grasp of the previous result is possible by realizing that the Kf B orbits and (|27|) are 
the flow lines respectively induced by the gradients of —W H and W 11 . Because 

W(^) = ( - 2 %) , hess"(^) = ( °) , ta"^) ' 1 " 



0^ and 0° are saddle points of W H whereas 4> B is a minimum. The flow lines of —W n accordingly start 
from cf> B and reach the cuspidal point (f> F at x = a, see Figure 9b. Then, new flow lines, now induced 
by W H , depart from F and go to ^ F is a conjugate point to B of this congruence of trajectories, 
which are thus unstable according to Morse Theory. 

We end this sub-Section by comment on an intriguing link with supersymmetry: our bosonic (1+1)- 
dimensional model admits two non-equivalent supersymmetric extensions, [2U- Acting on static on shell 
superfields of the form 

1 - dW 

$ a [x, e] = <p a (x) + er{x) - -ee— 

2 d(p a 

one finds two independent sets of M = 1 supercharges: 

nU - a d<t>a 4. ,/>« ft dwI n n rh -ft d<t)a + ,/,« ft dwU 

OU -ft d ^+ lh a n dWl Q U, -ft d ^+ ,,,a a ^ 

I are Grassman Majorana spinors that span the odd part of superspace; ifj a = I , \ ) 
t> 2 J \ 1P2 J 

are the two Majorana spinor fields, superpartners of the bosonic fields. Double solitary waves are classical 

BPS states, henceforth stable, because they are annihilated by a combination of the Q 1 supercharges: 

Q{^a = (Q{ + Qi)^a = (01 + d ^ + W + ^2) ~ (01 + 02)^~ = 

if, moreover, iftf+ip?; = 0. Quadruple solitary waves, however, are non BPS because Q^4> a + for them. 
As a curious fact, despite of having two sets of M = 1 supercharges this model does not admit J\f = 2 
supersymmetry because W(4>, <j>) = W^cfta) + iW n (4> a ) is non holomorphic. 



4 Adiabatic motion of non- linear waves 

In sum, the Kg (a, b) and K.f B (a, b) solitary wave families depend on a £ R, a translational parameter 
setting the center of the kink, and b £ R, the parameter measuring the distortion of the kink shape from 
a single lump. The space of kink solutions is thus the (a, b) 6 IR 2 -plane, with the basic kinks living in 
the boundary circle cM 2 ~ at infinity. In fact, the moduli space of double solitary waves is only the 
upper half-plane because invariance under 4>2 — ► —4>2 comes from invariance under b — ► — b for the 
kinks. 

We now ask the following question: Are there solitary wave solutions to the full wave equations 
©-® such that changes in shape take place? For normal kinks, the dynamics is merely dictated by 
Lorentz invariance and thus characterized by shape invariance. We shall address the difficult issue of the 
evolution of composite solitary waves within the framework of Manton's adiabatic principle, see [22] and 
geodesies in the moduli space determine the slow motion of topological defects. In this scheme, the 
dynamics arises from the hypothesis that only the parameters of the moduli space a and b depend on 
time. Moreover, the kink evolution (f> K (x; a(t), b(t)) is slow enough to guarantee that the static differential 
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equations will be satisfied for every given t in a good approximation. Therefore, the field theoretical 
action JQ) reduces to 



G 



J 1=1 



IK PtJK 



at at 



dt 



1 da da da db 1 db db 

2 y v ' y dtcft ycm ' ; tfttft 2 y v ' ; eft eft 



where 



9aa(a,b) = / 



tyf 9cg 
9a 9a 



9bb{x,b) = ^ 



dx- 



96 96 



i=i " - ,JU 

will be understood as the components of a metric tensor. 

We think of S G as the action for geodesic motion in the kink space with a metric inherited from the 
dynamics of the zero modes. The geodesic equations coming from the action S G determine the behaviour 
of the variables a(t) and b(t) and therefore describe the evolution of the lumps 4> K (x;a(t),b(t)) within 
the configurations of a specified kink family. 



4.1 Adiabatic evolution of double solitary waves 



The double solitary waves K^°(a,c)) can be seen as two particles moving on a line. Geodesic motion in 
the (a, 6)-plane would correspond to the slow-speed dynamics of the the two-particle system described in 
terms of the center of mass and relative coordinates a and b. 

For the value a = the formula ()24|) describing the K^ kinks is simpler and the tensor metric 
can be written explicitly. By changing variables to z = e x , the integrals become of rational type, and we 
obtain: 



g aa (a,b) = ^ 

9bb(a,b) 



g ab (a, b) = -gbb(a, b) = g ba (a, b) 



(1 - 



arccot 



b 2 ^T^j 



(35) 



From the graphics in Figure 10 we see that the metric induced in the (a, o)-plane by the zero mode 
dynamics has good properties: g aa > 0, g bb > 0, and, g = g aa g bb - g 2 b > 0. 



1 .l^Shbl b) 



0.5k?( b) 




3-2-1 1 2 £ 3 



Figure 10: Graphics from left to right of g ab , g b b and g = g a ag b b — g 2 ab as functions of b. 

In fact, metrics of the general form 

ds 2 = ^da 2 + 2F(b)dadb + G(b)db 2 

with both G(b) and g(b) = \G(b)-F 2 {b) greater than zero for all b, are always flat. Only two Christoffel 
symbols are non-null: 



r b _ 1 ( RF(h ) dF dG 
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Accordingly, the only relevant component of the curvature tensor is zero: -R^ab = 0. It is therefore 
possible to write the metric in the form 

ds 2 = da 2 + db 2 

by means of the isometric transformation (a(a, b), b(a, b)) if the new coordinates satisfy the PDE system: 



da\ ( db\ 1 da da db db (da 
~d~a) + \da~) = 4 ' d^"db + d^'db = ^' ' \db 



I -<*> 



(36) 



Zero curvature ensures the integrability of (|3f)|) and solutions for (a, b) can be found. In particular, 
restriction to the form (a = + a(6), 6 = ^= + 6(6)) provides the solution: 

a = V2 f db' [F(b') ± yW)] , b = V2 [ b db' [F(b') =F T^fe 7 )] 
jo ■/ o 



Geodesies are thus straight lines in the (a, b) variables: 

a(t) = ait + a 2 , 
where 01,02,61,62 are integration constants. 
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Figure 11: Evolution of the moduli variables and evolution of the lumps in the sector AO. 

Back in the original coordinates, a typical (a(t),b(t)) geodesic is depicted in Figure 11 (left). The 
b(t) curve shows that the relative distance between the two basic kinks decreases with time until they 
become fused -b = 0- at t = 50. Then, the distance starts to grow towards a complete splitting of 
the Jif B and Kf kinks in the remote future. The evolution of the center of mass is unveiled by the 
a(t) curve: it remains at x = while the two basic kinks are far apart. Shortly before the collision, a 
starts to grow and thereafter continuous to increase. The two kinks rebound, the center of mass being 
closer to the heaviest kink. The picture is clearer in the plot of the energy density evolution along this 
geodesic shown in Figure 11 (right). Shortly after t = —00, the geodesic starts from a point in the moduli 
space corresponding to one K± B kink and one Kf° kink very far apart. The two basic lumps begin to 
approach each other as \b\ decreases, distorting their shapes when they start to collide. Later, a single 
K^°(0) lump is formed when b(t) reaches (see Figures 5 and 11 (right)). At this point, the two lumps 
bounce back following the above-described motion in reverse, i.e., the separation between the two lumps 
increases, running asymptotically towards the configuration of two infinitely separated kinks in the C±° 
sector. Thinking of the graphics in Figure 5 as the stills of a movie, the adiabatic dynamics of the double 
solitary waves prescribes the speed at which the film is played. 
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4.2 Meta-stable slow motion of quadruple solitary waves 

We now perform the same analysis in the C BB sectors, again for the value a = The metric 

1 



h aa (a,b) = 2hbb(a,b) 



h ab (a,b) 



1 

32 



10 + 



(e fc -l)^ 



7rSign(e b — 1) — 2arctan 



1 



Ve b - 1 



(l-e 6 )^ 



-7rSign(l — e ) + 2arctan 



Vl - e b 



is of the general form studied in the previous sub-Section. 
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Figure 12: Graphics from left to right of h a b and h = h aa hbb — h 2 ab as functions of b. 



The geodesies are easily found and plotted using Mathematica in Figure 13 (left). The evolution of 
the energy density along this geodesic is also shown in Figure 13 (right). The motion begins at t = — oo 
with three lumps very far away from each other along the real line , arranged Kf°, K^°, Kf B 
configuration, ordered by increasing x. Then, the lumps start to approach the kink until they 

coincide, forming the K^ B (0) configuration. Later, the three lumps start to split in such a way that the 
Kf° kink overtakes the and Kf B lumps whereas the latter kink is also passed by the kink. 
This motion naturally proceeds towards Kf~ B , K^°, Kf° configurations with growing distances between 
the basic lumps. Note that there are no problems in this sector at b = ±oo; contrarily to the g- metric 
the /i-metric does not vanish at infinity. This means that some lumps could disappear smoothly as an 
effect of this low-energy scattering. 



-5 





b (t). 


t 


20 


P*- fin 


fin i nn 
a(* ) 




Figure 13: Evolution of the moduli variables and evolution of the lumps in the sector C± B . 

In fact, in Section §3 we argued that the K BB solitary waves are unstable by demonstrating the 
existence of a negative eigenfunction of the Hessian. Fluctuations in this direction of the configuration 
space will destroy the gentle evolution between the solitary waves described above. Closer understanding 
of the non-adiabatic process triggered by the negative fluctuation is possible by starting from a K BB (b) 
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lump with large b: there are one K± B and one Kf° lumps, respectively near x = — oo and x = +00, 
together with a kink sitting at the origin x = 0. The zero mode 



ycosh(x 



of the Hessian at K^°(0) provides a very good approximation to the negative eigenfunction of the Hessian 
at K BB (b — * 00) because this latter solitary wave tends asymptotically to K^°(0), ip^ has no nodes, 
and the topological barriers in the C±° and C BB sectors are different. Moreover, the ipQ fluctuation 
causes the splitting of the lump and the evolution reveals the motion of four lumps: two and 
two K BO kinks. Of course, this configuration does not belong to the static kink family 1)28(1 : the negative 
fluctuation suppresses the overlapping of two basic lumps. 

This picture is confirmed by solving via numerical analysis the Cauchy problem for the PDE system 
©-® with initial conditions 

<P(x,0) = Kf B (30)(x) , ^(x,0) =e^\x) 

and Dirichlet boundary conditions, see Figure 14. It seems that the two Kf° attract each other but 
there is repulsion between the two K± B kinks. 




Figure 14: Evolution of a ~K.4 BB (b) solitary wave after tfi^ fluctuations. 



4.3 Perturbing the system: non-linear wave bound state 

Finally, we add an small perturbation to our system. The basic kinks Kf and move adiabatically 
along the geodesies of the metric induced by the kinetic energy on the kink moduli space in the C±° 
sector. We choose to induce forces between the basic kinks by the following modification of V(cj>i, fa)' 

V'(fa,fa) = V(fa,fa) + ev(fa,fa) , v{fa,fa) = ( / ) lUl + j^-l 

Here, e is an small parameter and our choice is such that the zeroes of V are also zeroes of V 1 ; thus, 
V also comply with Coleman's theorem and enters into the class of admissible deformations of Vcsh in 
(l+l)-dimensions. Moreover, the Kf B and Kf° kinks survive this perturbation without any changes: 
they are still solitary wave solutions with the same energy. The solitary wave, however, needs to be 
adjusted: 

4>{x, a) = (-l) a ^-=\]l + (-l^tanhfyT + 2e(x - a)] , a,/3 = 0,l 




18 



is a solitary wave solution of the perturbed system with higher energy: £'{¥^2 AO ) = + 2e- The 
perturbation is chosen in such a way that the kink energy sum rule (|25|) is broken. 

All the other solitary wave are destroyed by the perturbation. Nevertheless, the adiabatic approxi- 
mation can still be used to describe the slow motion of the basic lumps. For small enough e, e « 1, the 
solitary waves (|24jl are approximate solutions of the perturbed system; therefore, we take the parameters 
a and b as the moduli variables, even though the perturbation forbids neutral equilibria in the moduli 
space. Plugging the approximate solutions into the action of the perturbed model we obtain: 



1 da da da db 1 db db 



dt 

The perturbation induces a potential energy, plotted in Figure 15 (left) 

/•oo 

u{b)= / dxv(K 2 AO (b)(x)) 



(37) 



and the motion ceases to be geodesic; u(b) is a potential barrier with maximum height reached at b = 
that goes to zero for b = ±oo. The equations of motion coming from (|37|) determine the evolution of the 
moduli variables a(t) and b(t) with time. 

In Figure 15 (center) the plot of both a and b as a function of time is shown for a solution of the 
equations of motion with e = 0.1 and initial conditions a(0) = 0, 6(0) = -2.004, a'(0) = 0, b'(0) = 0.001. 
The surprising outcome is that b seems to be quasi-periodic in time; this is a Sisyphean movement, 
crossing over and over again the top of the hill! 4 

In Figure 15 (right) the evolution of the energy density along this quasi-periodic orbit is drawn. We 
observe that the net effect of the perturbation plus the non-trivial metric is the upsurge of an attractive 
force between the two basic lumps Kf° and Kf B , giving rise to a bound state similar to the breather 
mode in the sine-Gordon model. 





0.3 


u(b) 


2 






b (t > 






.25 
















.1 
. 15 




1 






a It) > 






U.l 










60 


y\/ io 




OF. 05 




-1 

b 










A 


-2 


2 4 -2 












Figure 15: Plot of u(b), evolution of the moduli variables and evolution of the lumps in the deformed 
model. 
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